Effects of Phase Transition induced density fluctuations on pulsar dynamics 
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Abstract 


We show that density fluctuations during phase transitions in pulsar cores may have non-trivial effects on pulsar 
timings, and may also possibly account for glitches and anti-glitches. These density fluctuations invariably lead 
irfo non-zero off-diagonal components of the moment of inertia, leading to transient wobbling of star. Thus, 
’“Accurate measurements of pulsar timing and intensity modulations (from wobbling) may be used to identify the 
^^pecihc pattern of density fluctuations, hence the particular phase transition, occurring inside the pulsar core, 
^^hanges in quadrupole moment from rapidly evolving density fluctuations during the transition, with very short 
pfime scales, may provide a new source for gravitational waves. 
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I. INTRODUCTION 


The core of an astrophysical compact dense object, snch as a nentron star, provides physical conditions 
where transition to exotic phases of qnantum chromo dynamics (QCD) [T] may be possible. Snperflnid 
phases of nucleons are also believed to exist inside neutron stars with vortex depinning associated with 
glitches (though it may not provide a viable explanation for anti-glitches). In this paper we propose a 
technique to probe the dynamical phenomena happening inside the neutron star, which may also account 
for glitches and anti glitches in a unihed framework. 

We consider density fluctuations which invariably arise in any phase transition. We show that even with 
relatively very small magnitudes, these density fluctuations may be observable with accurate measurement 
of pulsar timings which can detect very minute changes in the moment of inertia (MI) of the pulsar. This 
provides a very sensitive probe for density changes and density fluctuations (especially due to formation 
of topological defects) during phase transitions in the pulsar core. Non-zero off-diagonal components 
of moment of inertia arising from density fluctuations imply that a spinning neutron star will develop 
wobble leading to modulation of the peak intensity of pulses (as the direction of the beam pointing 
towards earth undergoes additional modulation). This is a unique, falsihable, prediction of our model, 
that rapid changes in pulsar timings should, most often, be associated with modulations in changes in 
peak pulse intensity. It is important to note that the vortex de-pinning model of glitches is not expected 
to lead to additional wobble as the change in rotation caused by de-pinning of vortex clusters remains 
along the rotation axis. Density fluctuations will also lead to development of rapidly changing quadrupole 
moment which can provide a new source for gravitational wave emission due to extremely short time scales 
involved (despite small magnitude of this new contribution to the quadrupole moment) . 

The effect of phase changes on the moment of inertia has been discussed in literature. For example, 
moment of inertia change arising from a phase change to high density QCD phase (such as to the QGP 
phase) is discussed in ref. [2]. In the scenario of ref. [2], the transition is driven by slow decrease in the 
rotation speed of the pulsar, leading to increasing central density causing the transition as the central 
density becomes supercritical. It is assumed that as the supercritical core grows in size (slowly, over the 
time scale of millions of years), it continuously converts to the high density QGP phase (even when the 
transition is of hrst order). Due to very large time scale, the changes in moment of inertia are not directly 
observable, but observations of changes in the braking index may be possible. 

Our work differs from these earlier works primarily in our focus on the rapidly evolving density fluctu¬ 
ations arising during the phase transitions. This has not been considered before as far as we are aware. 
(Though effects of density inhomogeneities in a rotating neutron star have been discussed [3].) Density 
fluctuations inevitably arise during phase transitions, e.g. during a hrst order transition in the form of 
nucleated bubbles, and may become very important in the critical regime during a continuous transition. 
The density huctuations arising from a phase transition become especially prominent if the transition 
leads to formation of topological defects. Extended topological defects can lead to strong density huc¬ 
tuations which can last for a relatively long time (compared to the phase transition time). It is obvious 
that such randomly arising density huctuations will ahect the moment of inertia of the star in important 
ways. Most importantly, it will lead to development of transient non-zero oh-diagonal components of the 
moment of inertia, as well as transient quadrupole moment. Both of these will disappear after the density 
huctuations decay away and the transition to a uniform new phase is complete. Net change in moment 
of inertia will have this transient part as well as the hnal value due to change to the new phase. It seems 
clear that this is precisely the pattern of a glitch or anti-glitch where rapid change in pulsar rotation is 
seen which slowly and only partially recovers to the original value. As we mentioned above, transient 
change in quadrupole moment will be important for gravitation wave emission. It is important to realize 
that the net change in MI (as discussed previously, e.g. in m) is only sensitive to the diherence in the 
free energies of the two phases, and cannot distinguish different types of phase transitions. In contrast. 
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density fluctuations arising during phase transitions crucially depend on the nature of phase transition, 
especially on the symmetry breaking pattern (e.g. via topological defects). Identification of these density 
fluctuations via pulsar timings (and gravitational waves) can pin down the specific transition occurring 
inside the pulsar core. 

We consider the possibility of rapid phase transition in a large core of a neutron star. The scenario 
of slow transition as discussed in j2| is applicable for slowly evolving star (e.g. by accretion), with a 
transition which is either a weak first order, or a second order (or a crossover). If the transition is 
strongly first order then strong supercooling can lead to extremely suppressed nucleation rate, with 
transition occurring suddenly after the supercritical core becomes macroscopic in size. (This is similar 
to the situation of low nucleation rate leading to nucleation distance scales of the order of meters in the 
early universe as discussed by Witten |1].) Rapid transition in a large core of neutron star is naturally 
expected during the early stages of evolution of neutron star due to rapid cooling of star. It could also 
be driven by rapid accretion on a neutron star, coupled with a first order transition. 

We will not discuss any detailed scenario of such a rapid transition in this paper, and rather just note 
the possibility of such transitions. Our focus will be that once such a rapid transition occurs, what are 
its observational implications. Clearly, change in phase will lead to net change in MI of star, affecting 
its rotation. This has been discussed earlier, e.g. in [2]. We will focus on additional effects associated 
with presence of density fluctuations which lead to qualitatively new effects, not included in the earlier 
investigations. These effects are, a transient component of change in MI, development of the off-diagonal 
components of MI, and a transient, rapidly evolving quadrupole moment. 


II. EFFECTS OF DENSITY FLUCTUATIONS DUE TO BUBBLE NUCLEATION 

A rough estimate of change in the moment of inertia due to phase change can be taken from ref. [2] 
using Newtonian approximation, and with the approximation of two density structure of the pulsar. If 
the density of the star changes from pi to a higher density p 2 inside a core of radius Rq, then the fractional 
change in the moment of inertia is of order, 

^ ~ _ ij:^ 

/ 3>i 

Here R is the radius of the star in the absence of the dense core. For a QCD phase transition, density 
changes can be of order one. We take the density change to be about 30% as an example. If we take the 
largest rapid fractional change in the moment of inertia of neutron stars, observed so far (from glitches), 
to be less than 10“®, then Eq.(l) implies that Ro < 0.3 Km (taking R to be 10 Km). For a superfluid 
transition, we may take change in density to be of order of superfluid condensation energy density ~ 
0.1 MeV//m^ (see, ref.[5]). In such a case, Rq may be as large as 5 Km. These constraints on Rq arise 
from observed data on glitches/anti-glitches. These estimates may also be taken as prediction of possible 
large fractional changes in the moment of inertia (hence pulsar spinning rate) of order few percent when a 
larger core undergoes rapid phase transition. For example, Rq may be of order 2-3 km for QCD transition 
(from estimates of high density core of neutron star [B]), or it may be only slightly smaller than R for 
superfluid transition. 

We now discuss density fluctuations during phase transitions. First we consider fluctuations arising 
simply from a first order transition. Consider nucleation of relatively large number of bubbles (several 
thousand) inside the supercritical core. This is possible due to nonuniformities, even of purely statistical 
origin (e.g. from fluctuations in temperature 0 ). 
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A. Parameters for bubble nucleation and Results 


We simulate random nucleation of bubbles filling up a core of size Rq (= 300 meters). Effects of bubble 
nucleations will be characterized in terms of the following parameters. 

Bubble radius tq will be taken to vary from 5 meters to 20 meters, with bubble separation being of 
same order as bubble size (close packing). 

Density change in bubble nucleation is taken to be about the nuclear saturation density ~ 160 
MeV//m^ (e.g. for QCD scale). 

We hnd that density fluctuations lead to fractional change in MI, A//J ~ 4 x 10“® for tq = 20 meters 
implying similar changes in pulsar timings. Change in MI remains of same order when ro is changed 
from 20 meter to 5 meter. Due to random nature of bubble nucleation, off-diagonal components of the 
moment of inertia, as well as the quadrupole moment become nonzero and the ratio of both to the initial 
moment of inertia are found to be of order Ry/Io — Q/Iq — 10“^^ — 10“^°. 

This aspect of our model is extremely important, arising entirely due to density fluctuations generated 
during the transition. As these density fluctuations homogenize, hnally leading to a uniform new phase 
of the core, both these components will dissipate away. The off-diagonal component of moment of inertia 
will necessarily lead to wobbling (on top of any present initially), which will get restored once the density 
fluctuations die away. This will lead to transient change not only in the pulse timing, but also in the pulse 
intensity (as the angle at which the beam points towards earth gets affected due to wobbling). We again 
emphasize that the conventional vortex de-pinning model of glitches is not expected to lead to additional 
wobble as the change in rotation caused by de-pinning of vortex clusters remains along the rotation axis. 
Thus, the presence of intensity modulations associated with a glitch can distinguish between our model 
and the vortex de-pinning model. 

Generation of quadrupole moment has obvious implication for gravitational wave generation. One 
may think that a quadrupole moment of order Q/Jq ~ 10“^° is too small for any significant gravity wave 
emission. However, note that the gravitation power depends on the (square of) third time derivative of 
the quadrupole moment [8]. The time scales will be extremely short here compared to the time scales 
considered in literature for the usual mechanisms of change in quadrupole moment of the neutron star. 
Here, phase transition dynamics will lead to changes in density fluctuations occurring in time scales of 
microseconds (or even shorter as we will see below in discussions of topological defects generated density 
fluctuations). This may more than compensate for the small amplitude of quadrupole moment and may 
lead to these density fluctuations as an important source of gravitational wave emission from neutron 
stars, as we will discuss below. 


III. DENSITY FLUCTUATIONS FROM TOPOLOGICAL DEFECTS 

Topological defects form during spontaneous symmetry breaking transitions via the so called Kibble 
mechanism |9]. These defects can be source of large density fluctuations depending on the relevant energy 
scales, and their formation and evolution shows universal characteristics, (e.g. scaling behavior). This 
may lead to reasonably model independent predictions for changes in MI, and quadrupole moment and 
subsequent relaxation. As bubbles, strings, domain walls, all generate different density fluctuation, with 
specihc evolution patterns, high precision measurements of pulsar timings and intensity modulations 
(from wobbling) and its relaxation may be used to identify different sources of fluctuations, thereby 
pinning down the specihc phase transition occurring. Specihc phase transitions expected inside pulsars 
lead to diherent types of topological defects. Important thing is that a random network of defects will 
arise in any phase transition, and resulting defect distribution can be determined entirely using the 
symmetry breaking pattern. For example, superhuid transition leads to formation of random network 
of vortices. Conhnement-deconhnement QCD transition can lead to formation of a network of domain 
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walls and global strings arising from the spontaneous breaking of Z(3) center symmetry [10]. QCD 
transition may also give rise to only global strings, e.g. in the color flavor locked (CFL) phase with 
SU{3)c X SU{3)l X SU{3)r x U{1)b symmetry (for 3 massless flavors) is broken down to the diagonal 
subgroup SU(3)c+l+r x Z 2 DP- As the resulting density fluctuations are mostly dominated by the nature 
of defect (domain walls, or strings, or both), we will model formation of these different defect networks 
in the following using correct energy scale and try to estimate resulting density fluctuations. 


A. Results of model simulation of defect network 

We model formation of U(l) strings and Z 2 domain walls by using correlation domain formation in a 
cubic lattice, with lattice spacing ^ representing the correlation length, as in ref. m- It is not possible 
to carry out these simulation covering length scales of km (for star) to fm (QCD scale). Hence these 
simulations are necessarily restricted to very small system sizes. Each lattice site is associated with an 
angle 6 randomly varying between 0 and 27i (to model U(l) global string formation), or two discrete 
values 0,1 (when modeling Z 2 domain wall formation). For string case, winding of 6 on each face of the 
cube is determined using the geodesic rule. For a non-zero winding, a string segment (of length equal 
to .^ ) is assumed to pass through that phase (normal to the phase). For domain wall case, any link 
connecting two neighboring sites differing in Z 2 value is assumed to be intersected by a planar domain 
wall (of area and normal to the link). The mass density (i.e. mass per unit length) of the string was 
taken as 3 GeV/fm, and the domain wall tension is taken to be 7 GeV/fm^. These values are taken as 
order of magnitude estimates from the numerical minimization results in ref. [lO] for the pure gauge case. 
(Note that logarithmic dependence of global strings on inter-string separation may lead to much larger 
density fluctuations than considered here.) 

We consider spherical system of size R and conhne defect network within a spherical core of radius 
Rc = ^R- This is in view of the constraints on the supercritical core size Rc being of order 0.3 km for a 
neutron star with radius i? = 10 km. Of course, in our simulations, R is extremely small, with maximum 
value of 4000 fm. For ^ ~ 10 fm, we hnd the resulting value of ^ ~ 10“^^ — 10“^^ implying similar 
changes in the rotational frequency. Here RR = 1,2,3 are the three diagonal values of the MI tensor. 
As we increase Rc from 5 to about 400 we hnd that the value of ^ stabilizes near 10“^^ — 10“^^ as 
shown in the table below. This change in ^ amounts to change in the number of string and wall segments 
by a factor of 10®. This gives a strong possibility that the same fractional change in the MI may also 
be possible when R is taken to have the realistic value of about 10 Km, especially when we account 
for statistical huctuations in the core. For the formation of domain walls we hnd fractional change in 
MI components (as well as quadrupole moments) to be larger by about a factor of 40. With accurate 
measurements, these small changes may be observable. Also, for a larger core size undergoing transition 
these changes can be larger. We are only presenting change in MI due to transient density huctuations 
during phase transition, which as we see, can have either sign. (As we discussed above, the net change in 
MI will include the very large contribution of order 10“® due to net phase change of the core [2].) This 
suggests that the phase transition dynamics may be able to account for both glitch and anti-glitch events 
(with associated wobbling from oh-diagonal components leading to pulse intensity modulation). 

We now consider superhuid transition. A rapid superhuid transition could occur after transient heating 
of star (either due to another transition releasing latent heat, or due to accretion etc). (In this work 
we neglect any possible ehect of star rotation on this mechanism of random vortex formation which 
is expected to lead to a much denser network than the one arising from star rotation.) We take the 
vortex energy per unit length to be 100 MeV/fm and correlation length for vortex formation of order 
10 fm (ref. [5|). Table I shows that the string induced transient fractional change in MI is of order 10“^® 
(compared to net fractional change in MI of order 10“® as discussed in Sect.H). Ratios of the quadrupole 
moment and off-diagonal components of MI to the net MI of the pulsar are also found to be of order 
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10 -io_ transient change in the MI decays away when the string system coarsens, thereby restoring 
less than few percent of the original valne. We note that this pattern (and nnmbers) are similar to that 
of a glitch (or anti-glitch). 


TABLE I: Fractional change of various moments of the pulsar caused by inhomogeneities due to defects, with the 
correlation length ^ = 10 fm. For QCD scale strings, the string tension is taken as 3 GeV/fm, while the QCD 
Z(3) wall tension is taken as 7 GeV/fm^ (from simulations in ref.|10]). For the superfluid vortices, the energy 
per unit length is taken to be 100 MeV/fm [5]. 
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B. Field theory simulations for QCD transition 

As a further support for these estimates, we have also carried out held theory simulations of 
conhnement-deconhnement (C-D) QCD transition using effective held theory Polyakov loop model. This 
leads to spontaneously broken Z(3) symmetry (for the SU(3) color group) in the QGP phase giving rise to 
topological domain wall defects in the QGP phase and also string defects forming at wall junctions. We 
carry out a held theory simulation for the C-D transition using a quench (quench is used for simplicity 
as only domain formation is relevant here), see ref.[Tn] for details. It is not possible to carry out held 
theory simulation covering length scales of km (for star) to fm (QCD scale). Hence these simulations 
are necessarily restricted to system sizes of tens of fm only. The physical size of the lattice is taken as 
(7.5 fm)^ and (15 fm)^. We use periodic boundary conditions and take a spherical region with radius 
Rc (representing star’s core) to study change of MI, with Rc = 0.4 x lattice size. We mention that the 
model of ref. [10] does not directly apply to the case of neutron star which has large baryonic chemical 
potential. However, the most relevant features of the model are formation of Z(3) defects with similar 
energy scale as in the neutron star case. Hence we use those simulations in to estimate defect induced 
density huctuations for the present neutron star case. We also use dissipation to relax density huctuations 
but total energy is kept hxed by adding the dissipated energy. Thus we only focus on re-distribution of 
energy in defect network and the background. For the net change in the MI, we will use the estimates of 
A/// from Sect.H (refs. [2]). QCD transition may also give rise to only global strings, e.g. in the color 
flavor locked (CFL) phase with SU{3)c x SU{3)l x SU{3)r x U{1)b symmetry (for 3 massless flavors) is 
broken down to the diagonal subgroup SU{3)c+l+r x Z 2 [T]. For this case, we modify the model studied 
in ref. [10] by removing terms which correspond to Z(3) structure of the vacuum manifold. This gives 
rise to string defects only without any domain walls, with energy scale of QCD. 

Plots in Fig.l show the time evolution of the resulting fractional change in the MI and the quadrupole 
moment of the core relative to the initial total MI for these held theory simulations. Here we have taken 
dense core of fractional size 0.3/10 (in view of transition in core size of 300 meter for a 10 km star). We 
add the MI of a shell outside the core so that total system size = jjh of the core size, and the shell has 
the same uniform density as the core. Fractional changes for all components of MI as well as the ratio 
of quadrupole moment to MI are found to be of similar order. We are only presenting change in MI due 
to transient density huctuations during phase transition, which as we see, can have either sign. (Net 
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FIG. 1: Fractional change in MI and quadrupole moment during phase transitions. (a),(b) correspond to lattice 
size (7.5 fm)^, and (c),(d) correspond to lattice size (15 fm)^ respectively. Plots in (a) and (c) correspond to the 
C-D phase transition with Z(3) walls and strings, while plots in (b) and (d) correspond to the transition with 
only string formation as for the CFL phase. 

fractional change in MI will include the very large contribution of order 10“^ due to net phase change of 
the core [ 2 ].) This suggests that the phase transition dynamics may be able to account for both glitch 
and anti-glitch events (with associated wobbling from off-diagonal components leading to pulse intensity 
modulation). As discussed above, rapid changes in the quadrupole moment will lead to gravitational 
wave emission. 


IV. GRAVITATIONAL WAVE GENERATION DUE TO DENSITY FLUCTUATIONS 


Despite the small values of quadrupole moments in Table I, the power emitted in gravitational waves 
may not be small due to very short time scales in the present case. The defect coarsening will be governed 
by microphysics with time scale being of order tens of fm/c. Even with extremely dissipative motion of 
strings, and much larger length scales, the change in quadrupole moment due to strings can happen in an 
extremely short time scale (during string formation, and/or during string decay), thereby boosting the 
rate of quadrupole moment change. Even a very conservative value of the time scale of microseconds (e.g. 
for nucleation of bubbles with nucleation sites few meters apart, as discussed above) for the evolution 
of density fluctuations will still be 1000 times smaller than fastest pulsar rotation time of milliseconds, 
resulting in huge enhancement in the gravitational power, as can be seen from the following expression 
for the power |8] 


1-3 




sec. 


Y 


( 2 ) 


Here Jq is the MI of the pulsar, and AQ is the change in the quadrupole moment occurring in time 
interval At. For the present case, we can take AQ/Iq of order — 10“^*^ from Table I (which is much 
smaller than the value of 10“®, typically used for deformed neutron stars). Even though AQ/Iq is very 
small here, the relevant time scale At is also very small. Thus, even with a conservative estimate of the 
time scale for the transition (e.g. bubble coalescence). At = 10“® — 10“® sec., the power in gravitational 
wave can be signihcant due to large enhancement from the factor. A shorter time scale may 

make this source as potentially very prominent for gravitational waves. For an estimate of the expected 
strain amplitude from a pulsar at a distance r, we use the following expression [8] 


h = 


Att^GAQP 

c^r 


( 3 ) 

With AQ/Iq of order 10“^° and a time scale for the transition (e.g. bubble coalescence). At = 
10“® — 10“^ sec. we can have h ~ 10“^^ — 10“^^ for a pulsar at 1 kpc distance. As we mentioned 
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above, the time scale for evolution of density fluctuations may be even shorter, leading to larger strain 
amplitudes (even if much smaller values of AQ/Iq are taken from Table I), and much larger power emitted 
in gravitational waves. Since the wave emission is only for a single burst, lasting for only duration At, 
net energy lost by the star remains small fraction of the star mass. 


V. CONCLUSIONS 

To summarize our results, we have shown that density fluctuations arising during a rapid phase 
transition lead to transient change in the MI of the star. Such density fluctuations in general lead to 
non-zero off-diagonal components of moment of inertia tensor which will cause the wobbling of pulsar, 
thereby modulating the peak intensity of the pulse. This is a distinguishing and falsihable signature of 
our model. The conventional vortex de-pinning model of glitches is not expected to lead to additional 
wobble as the change in rotation caused by de-pinning of vortex clusters remains along the rotation 
axis. We hnd that moment of inertia can increase or decrease, which gives the possibility of accounting 
for the phenomenon of glitches and anti-glitches in a unihed framework. Development of nonzero value 
of quadrupole moment (on a very short time scale) gives the possibility of gravitational radiation from 
the star whose core is undergoing a phase transition. Net change in MI (as discussed previously, e.g. 
in 0) is only sensitive to the difference in the free energies of the two phases, and cannot distinguish 
different types of phase transitions. In contrast, density fluctuations arising during phase transitions 
crucially depend on the nature of phase transition, especially on the symmetry breaking pattern (e.g. 
via topological defects). Identihcation of these density fluctuations via pulsar timings (and gravitational 
waves) can pin down the specihc transition occurring inside the pulsar core. This is an entirely new way 
of probing phase transitions occurring inside the core of a neutron star. Though our estimates suffer 
from the uncertainties of huge extrapolation involved from the core sizes we are able to simulate, to the 
realistic sizes, they strongly indicate that expected changes in moment of inertia etc. may be well within 
the range of observations, and in fact may be able to even account for the phenomena of glitches and 
anti-glitches. 

One aspect of glitches which may raise concern in our model is multiple occurrences of glitches. For 
vortex de-pinning model multiple glitches seem natural due to formation of vortices due to star rotation 
with glitches occurring with de-pinning of clusters of vortices. In our model, multiple occurrence of glitches 
will require multiple phase transitions. This is not very improbable, though, when the transition happens 
due to accretion from companion star causing heating, with subsequent cooling down. Essentially, in this 
case, the matter in the neutron star core will continue to remain in the vicinity of the phase boundary 
in the T — /i plane (in the QCD phase diagram). Accretion will move the system across the boundary 
by increasing /i, causing a phase transition, and subsequent cooling will move the system back through 
the phase boundary with another transition. Such processes could in principle repeat for certain neutron 
stars leading to multiple glitches. One should also allow the possibility that glitches could occur due to 
multiple reasons. Some glitches/anti-glitches could occur due to the model proposed here, that is due to 
phase transition induced density fluctuations, while other glitches could occur due to the conventional 
de-pinning of vortex clusters. Main point of our work is to emphasize that if and when a transition 
happens in the core of a neutron star, it invariably leads to density fluctuations which manifest itself in 
glitch/anti-glitch like behavior, along with other implications such as wobbling of star, gravitational wave 
emission etc. 

It is thus important to investigate this interesting possibility further. With much larger simulations, 
and accounting for statistical fluctuations of temperature, chemical potential etc. in the core, more 
dehnitive patterns of changes in moment of inertia tensor/quadrupole moment etc. may emerge which may 
carry unique signatures of specihc phase transitions involved. (For example, continuous transitions will 
lead to critical density huctuations, and topological defects will induce characteristic density huctuations 


depending on the specific symmetry breaking pattern.) If that happens then this method can provide a 
rich observational method of probing the physics of strongly interacting matter in the naturally occurring 
laboratory, that is interiors of neutron star. It will be interesting to see if any other astrophysical body, 
such as white dwarf, can also be probed in a similar manner. 
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